(An)Isotropic models in scalar and scalar-tensor cosmologies 
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We study how the constants G and A may vary in different theoretical models (general relativity 
with a perfect fluid, scalar cosmological models ("quintessence") with and without interacting scalar 
and matter fields and a scalar-tensor model with a dynamical A) in order to explain some observa- 
tional results. We apply the program outlined in section II to study three different geometries which 
generalize the FRW ones, which are Bianchi V, Vllo and IX, under the self-similarity hypothesis. 
We put special emphasis on calculating exact power-law solutions which allow us to compare the 
different models. In all the studied cases we arrive to the conclusion that the solutions are isotropic 
and noninflationary while the cosmological constant behaves as a positive decreasing time function 
(in agreement with the current observations) and the gravitational constant behaves as a growing 
time function. 



I. INTRODUCTION 

In a series of recent papers we have studied 

and compared how the "constants" G and A may vary 
in different theoretical frameworks for several metrics. 
These different theoretical frameworks are: general rel- 
ativity with time varying constants (TCV), scalar cos- 
mological models with non-interacting scalar and matter 
fields and TVC, and the last of the studied models is the 
usual scalar-tensor theory with a dynamical cosmological 
constant which seems to be the most natural theoretical 
model to study the possible variation of the gravitational 
and the cosmological constants. In those recent works we 
have been able to state and prove general results valid 
for all the geometries (all the Bianchi types as well as for 
the FRW geometries) within the context of self-similar 
solutions (SSS). We have focused our attention on this 
class of solutions since, as has been pointed out by Coley 
(0)' ttis self-similar models play an important role in de- 
scribing the asymptotic dynamics of the Bianchi models. 
A large class of orthogonal spatially homogeneous mod- 
els (including all class B models) are asymptotically self- 
similar at the initial singularity and are approximated by 
exact perfect fluid or vacuum self-similar power-law mod- 
els. In the same way, exact self-similar power-law models 
can also approximate general Bianchi models at interme- 
diate stages of their evolution and which is also impor- 
tant, self-similar solutions can describe the behaviour of 
Bianchi models at late times. Working under the hypoth- 
esis of self-similarity allows us to find exact power-law so- 
lutions in such a way that we may compare the obtained 
solution for each studied model. 

In this paper we extend our program by considering 
the important case of the interacting scalar and matter 
fields within the framework of scalar cosmological mod- 
els, and of course, we formulate the corresponding model 
with TVC. Therefore the aim of the present work consists 
in studying and comparing, by calculating the exact so- 
lution, how the constants G and A may vary in different 
theoretical models and with different geometries under 
the self-similar hypothesis. We apply the outlined pro- 



gram to study three different geometries which generalize 
the FRW ones, which are Bianchi V, VIIq and IX, since 
as pointed out in contrary to what is widely be- 

lieved, an almost isotropic cosmic microwave background 
(CMB) temperature does not imply that the universe is 
"close to a Friedmann-Lemaitre universe" and for this 
reason it is important to study anisotropic models. 

For this purpose, in section II we start by summarizing 
the main results that we employ in the next sections. In 
section IIA, we begin by setting the notation and explain- 
ing how to calculate a homothetic vector field (HVF). 
We also give the definitions of the some geometrical and 
kinematical quantities, such as the anisotropic and the 
deceleration parameters. In section IIB we outline the 
so called classical models. These models are the general 
relativity one with a perfect fluid and its TVC model de- 
rived from it. We show the general form of the solutions 
in each case. Section IIC is devoted to the scalar cos- 
mological models. We review the results stated in (Q) 
showing the admitted form for the scalar field and its 
potential for every homothetic solution (independently 
of the geometrical model). We study the following cases: 

• The standard scalar cosmological models with G' = 
0, i.e. where the gravitational constant is consid- 
ered as a true constant and its corresponding gen- 
eralization, i.e. a scalar model with G' ^ 0. 

• The non-interacting case of a scalar and matter 
fields and its TVC generalization. 

• As a part of our enlarging program, we now study 
the interesting case of an interacting scalar and 
matter fields. We compare three "a priori" (in 
principle) different approaches. The first one pro- 
posed by Maia et al (16]) where the authors do 
not split the resulting conservation equation follow- 
ing a thermodynamical approach. The second one 
has been proposed by Wetterich ([26i|). In this case 
the author consider the splitting conservation equa- 
tions and takes as the possible coupling between 
the scalar and matter field the function g"^ = Spm, 
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where the constant S must be negative. The third 
of the studied cases is the one proposed by BiUard 
et al where in analogy with the model pro- 
posed by Wetterich but following a different point 
of view the authors consider as a coupling function 
Q = 6Hpm, where, in this case the constant S must 
be positive in order to satisfy the second law of 
thermodynamics. 

• We conclude the study of the scalar cosmological 
models by considering the case of an interacting 
scalar and matter field with TVC. For this purpose 
we split the generalized conservation equation and 
consider as a possible coupling function Q = SHp,yi, 
with S > 0. 

In section IID we consider a general scalar-tensor the- 
ory with a dynamical cosmological constant, A. 

Once we have presented the program, then we ap- 
ply it, step by step, to three difi^erent geometries which 
are Bianchi V, Bianchi VIIq and Bianchi IX, since they 
may be considered as an anisotropic generalization of the 
FRW geometries. In section III we study the Bianchi V 
model, while in sections IV and V are studied the mod- 
els Bianchi VIIq and IX, respectively. We begin each 
section by deducing the metric through their Killing vec- 
tor fields (KVF) and outlining the Einstein tensor. Then 
we go next to study the classical, scalar and scalar-tensor 
models as we have described previously. In each case we 
calculate the standard solution, i.e. with the constants 
acting as true constants and its TVC solution in order to 
show how such hypotheses modify the standard solution. 
We end in section VI with a brief conclusions. 

In the appendix A we have included the proofs of the 
results stated in section II concerning the interacting 
scalar and matter models. We give a rigorous proofs of 
the results by using the Lie group method. We also ex- 
plore a generalization which allows us to consider differ- 
ent coupling functions between the scalar and the matter 
fields. In appendix B we prove the results for the case of 
an interacting scalar and matter fields with TVC. 



II. THE MODELS 

A. Self-similar solutions 

Throughout the paper, Ai will denote the usual 
smooth (connected, Hausdorff, 4-dimensional) spacetime 
manifold with smooth Lorentz metric g of signature 
(— , (see for example [2l|). Thus M is para- 
compact. A semi-colon and the symbol C denote the 
covariant and Lie derivative, respectively. We shall use 
a system of units where c — 1. For a metric, g, and for 
a vector field H G X{M), {'H_—'Hi{t,x,y,z)dxi)i^i the 
homothetic equation reads ( 23] and Igi]): 



In the case of the Bianchi models, this equation brings 
us to obtain that the scale factors behave as follows 

a = ao{t + tof' , b = bo{t + to)"' , d = do {t + tof' , 

(2) 

where {ai)^^^ G M^. In that follows we define, H — 
h{t + to) ^ , with h = ai + a2 + 03. In each of the stud- 
ied cases we will get restrictions on the constants (oj) .^-^ . 
We define the deceleration parameter as 

3 



(3) 



and the anisotropic quantities, A ([1|-[1|), and fj25j): 



■^=Jp^ (4) 

that give us a measure of the anisotropy. is the shear 
scalar and I3 = C°'^'^'^Cabcd, is the Weyl scalar. 



(1) 



B. The classical models 

The field equations (FE) read 

R^J - \R9^3 = ^^GT,j - Kg,j , =0, (5) 

where we consider vacuum solutions by making Tij = 
A = 0, and the classical perfect fluid solutions if A = 0. 
For a perfect fluid (PF) the energy-momentum tensor is 
deflned by: 

^" ^{P + P) UiUj + pg,j , (6) 

where p is the energy density of the fluid, p the pressure 
and they are related by the equation of state p = up, 
(w G (— 1, 1]) , Ui is the 4— velocity. From the conserva- 
tion equation, T^.j = 0, and taking into account that the 
scale factors are given by Eq. ^ , then it is easy to arrive 
to the conclusion that the energy density behaves as 

p = po(t + to)-2. (7) 

In order to take into account the variations of G and A 
we use the Bianchi identities (see for example [l5|-[l^|) 

= (S^GTy -Ag,,)'^ (8) 
which read: 

St:G[p' + p{l+uj)H] = -A' -SnG'p, (9) 

(' means time derivative) in our case we obtain (assuming 
the additional condition, T-.^ = 0) : 

p' + p{l+Lu)H = 0, and A' = -8nG'p. (10) 

It could be proven that the solution for this model (and 
for every Bianchi model) has the following form 

p = Poit + toy" , G^Go{t + toT-\ 
A = Aoit + toy\ (11) 
where a = /i (w + 1) G M, and po- Go, Aq G M. 



3 



C. Scalar models 

We consider the following cases: 

1.- For a scalar field (/>, the stress-energy tensor may 
be written in the following form (3): 



(12) 



where the energy density and the pressure of the fluid 
due a scalar field are given by 

P^^\^'^ + V{c^), = (13) 

and conservation equation now reads (the Klein-Gordon 
equation (KG)) 



h" + H(j)' + —V = 0. 



(14) 



In a previous work ( [21 ) we have proven that the unique 
form for the scalar field, 0, and the potential, V{(j)), com- 
patible with the self-similar solution (SSS) are given by 



(j) = ±^/a\-a{t + to), F = /3exp(=F 



(15) 



where a e M+, while ^ e M. 

2.- We would like to study how the gravitational con- 
stant varies when we are considering a scalar field. For 
this purpose, by using the Bianchi identity (G(t)T'^)^., = 
0, we get 



(16) 



which is the modified KG equation. In the G— var frame- 
work, in (0) we have proven that the main quantities 
behave as follows 

(/. = 0o(t + io)"", y-/3(t + to)"'^"+'\ 
G-Go(t + io)', (17) 

where a, 0Oi Go, /3 € M, and g — 2a. Note that they must 



verify Gi/)'' 



r 



GV « t . For a similar approach in 



the context of a FRW model with holographic dark en- 
e rgy with varying gravitational constant see for example 

Scalar and matter fields models. In this case the stress- 
energy tensor may be defined by, T = T"^+T'^, where T™ 
is defined by Eq. © and T"^ by Eq. (HH). We may con- 
sider two cases: The non-interacting and the interacting 
one (with their TVC versions). 

3.- In the non- interacting case we already know how 
each quantity must behave since we consider the conser- 
vation equations (divT™ = = divT"^) so they behave 



(18) 



= ±v^ln(i + io), V = l3exp[T 



This class of solutions are known in the literature as 
scaling cosmological solutions (Q and the references 
therein). 

4.- For the non-interacting case with G— var we con- 
sider, by taking into account the Bianchi identity, the 
following modified KG equation (see for details) 



Gp' + Gip + p)H ^ -G'p, 



(19) 



where p — p„i + p^f,, and p = Pm +P,p, which we may split 
into 

G' 

P'm + (Pn, + Pm) = --^Pm, (20) 

where we already know how each quantity behaves, 
namely 

<t>^4>o{t + toy", ^/ = /3(^ + ^o)"'^"+'^ 

G^Goit + tof", p = po(t-fio)-2(a+i)^ (22) 

5.- In the case of an interacting perfect fluid with a 
scalar field, interacting quintessence, we may consider 
several possibilities. In these models, it is considered that 
the scalar field decays into the perfect fluid and therefore 
it may be interpreted as a model with dark matter (DM) 
coupled to dark energy (DE). In order to analyze this 
case we consider three different approaches: 

1. Approach 1. ([l^) In this case the conservation 
equation reads 

p'„, + {u; + l)p,nH = -<!>' (acl>+^y (23) 

In the appendix we will prove that the behaviour 
of the quantities are given by Eq. In f[i6j) 

the authors interpret the scheme from a thermody- 
namical point of view, a la Prigogine ([13]) i.e. a 
model with matter creation. 

2. Approach 2. ([26]) In this case the conservation 
equation is split as follows 



Pni + (w + 1) PmH 
dV 
~d4 



□0- 



(24) 
(25) 



where q'^ = Sp,n- Therefore we have obtained, see 
appendix, that the quantities behave as Eq. ()18|) . 
Setting S — 0, we regain the non-interacting case. 
The only constrain for the model is (5 < 0. The 
author considers the possible couplin g of the scalar 
fleld to matter from previous works (p7|). 

3. Approach 3. (4]) In this case the conserva- 
tion equations ((T™)^.^. = (5i?p„, and (T"^)^.^. = 
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-SHpra,) read 

p'„^ + + 1) p,nH = Q, 



0-^+ -7T 



-Q, 



(26) 
(27) 



where Q = SHpm- Again, in the appendix we will 
show that the behaviour of the quantities are given 
by Eq. (fT8|) . As above, the energy of the scalar 
field is transferred to the matter field. The only 
restriction is to assume (5 > 0, in such a way that 
the second law of thermodynamics is verified f [l8|). 
Note that if (S = 0, then the non-interacting case is 
regained. The authors obtain the interacting term, 
6Hpm, by taking into account the conformal trans- 
formation that relates the scalar-tensor theory to 
the scalar one. 

In the appendix we will discuss the most general form 
of the function Q, in the framework of the self-similar 
solutions. 

6.- To end, we consider the case of interacting fluids 
within the framework of G— varying. For this model the 
modified KG equation reads 

p'^ + + ^0' + H {^'^ + {pra + P,n)) = 



G' 

G 1 P" 



V 



where we may split (as above) it in order to get 

G' 

P'm + (Pni + Prn) H + — p„i = Q, 
dV\ G' 



with for example 



Q = SHpr, 



(28) 

(29) 
(30) 

(31) 



In this case, the solution takes the form (see the ap- 
pendix for details): 

= 00 (t + to)-" , v{t) = p{t + toy'^"+'^ , 

G = Go (i + io)'" , Prn^Poit + io)"'^"+'^ • (32) 
For an alternative point of view see for example ([H). 



D. Cosmological models with dynamical A in 
scalar-tensor theories 



where g = det(gy), G* is Newton's constant, Sng is 
the action for the nongravitational matter. The arbi- 
trary functions uj (0) and A (0) distinguish the different 
scalar-tensor theories of gravitation, A ((f)) is a potential 
function and plays the role of a cosmological constant, 
and UJ (d)) is the coupling function of the particular the- 
ory ([281). 

The explicit field equations are 



87r„ . , , s UJ 
—T,j + A {4>) g,, + — 



gijOq 



(34) 



(3-I-2W ((/)))□(/. = 87rr- 



duj , , , ^ / dA . , , , 
— (/..,0^'-20 ^—~Ai(f>) 



(35) 

where, T — T/, is the trace of the stress-energy tensor, 
defined by Eq. (jH]). The gravitational coupling Gce{t) is 
given by 



Gesit) 



/ 2cj + 4 \ G, 
\2uj + 3) <j){ty 



(36) 



In a recent work (see [1] ) we have proven that the self- 
similar solution admitted for the EE p4ll35l) have the 
following form 



= 0o(t + to)", A((^)=Ao0' 



(37) 



with fi + a ^ 2, therefore A{t) = Ao{t + to) . The 
Brans-Dicke parameter is constant, w (0) — const., and 
p = po{t + to)~" , a = (1 + j) h. Note that in this case 
we have changed the notation and now the parameter of 
the EoS is 7, so p — jp. 



III. BIANCHI V MODELS 

A. The metric 

We consider (see for example ) the following Killing 
vectors fields (KVF) 

6 = 9y, ^2 = dz, ^3 = (9a. -I- mydy + mzdz, (38) 



We start with the action for the a general scalar-tensor 
theory of gravitation with A {(jj) , 



1 



167rG* 
Sng, 



-9 



<^ii-^^^M4M.+2</>A(0) 

(33) 



where m e M\ {0} , such that, — C^^^k, i-e. 

[6,6] = 0, [6,6] = Ki, and [6,6] = Ki- With 
these KVF we obtain the following metric 

ds^ = -dt^ + a{t)^dx^ + b{t)'^e~'^"'''dy'^ + d{t)^ e-'^"''' dz^ . 

(39) 
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The Einstein tensor reads 



b'a' d'a' d'b' 3to2 ^ 

ba da db a? 

a b d 

b" d" d'b' n?^ _ ^1 

b d db ^' 

d" a" d'a' toV _ 

~T ^ ~i 2~ ~ ""2 J 

d a da 

b a ba 



- r3 



(40) 
(41) 
(42) 
(43) 
(44) 



From Eq. ([T|) we obtain the following homothetic vec- 
tor field 

n^tdt + {l~ aa) ydy + (1 - as) zdy, (45) 

with the following constrains: 

a = aoit + to), b = boit + to)"" , d ^ do (t + to)"' , 

(46) 

where ai = 1 and (ai)j^2 ^ -1^' ^'^'^ therefore, H — 
h{t + to)~^ , with h = 1 + a2 + as. Note that we have 
obtained a non-singular solution. 



B. Classical solutions 

1. Vacuum solution 

The only possible solution is, 02 = 03 = 1, and m = 
±1, hence the solution is 

a = ao{t + to), b^bo{t + to), d^ do {t + to) , (47) 

with m — ±1. Therefore q — 0. With these results the 
metric yields 

ds^ = -dt^ + {t + to)^ {dx^ + 6=^2"=" + dz"^)) , (48) 

and therefore it admits more KVF, ^4 = —mzdy + mydz, 
^5 = -2mzd,^ - 2m^yzdy + {m^y^ - m'^z'^ + e^™^) 9^, 
and ^6 = 2mydx + {m?y'^ — m'^z^ — e^™^) dy + 2m?yzdz. 
This metric is the Milne form of a flat space-time (see 
Chapter 9, Eq. (9.8) of iM])- 



and Vm G (— 1,1)\{0}. With these results the metric 
collapses to Eq. (j48|) . it does not inflate, (7 = 0, and it 
is only valid for the EoS, a; = — ^, which is no strange 
since for example the SSS for the Bianchi I model is only 
valid for w = 1. A simple calculation shows us that the 
solution is isotropic since the anisotropic parameters A 
and vanish. From the DS point of view, the solution 
is stable ([131) since it is the FRW model with negative 
curvature. 



3. Time varying constants model 

In this case the behaviour of the solution is given by 
Eq. pn . The EE for this model are described by Eqs. 
(I40II44P with the conservation equations (ITU|) , where a = 
[lo + l)h, and h = (1 + a2 -f 03) , and the stress-energy 
tensor is defined by Eq. ([6]). Therefore we have found 
the next solution 



0,1 — a2 — as ~ 1, q ~ 0, a = 3 (a; -|- 1) . 



(50) 



(1 



(1- 



Go = ^ — -, ^, Aq = 3 (1 - m^) - \ 

with the following behaviour for G and A 



G; 



Ao 



decreasing e (—1, —1/3) 
constant if cj = — 1/3 
growing Vw G (—1/3, 1] 

negative Vw G (—1, —1/3) 
vanish if w = — 1/3 , 
positive Vw e (—1/3, 1] 



this solution is valid Vw e (-1, 1], m e (-1, 1) \ {0} . As 
above, with these results the metric collapses to Eq. pS)) 
and therefore it does not inflate. It is isotropic and valid 
Vw € (—1, 1]. As it is observed, if we set oj = —1/3, then 
the solution collapses to the above one with G = const. 
and A = 0. With regard to the behaviour of G and A, we 
only may say that "if we take into account the current 
observations ([l^) which suggest us that Aq > 0, then 
the solution is only valid Vw € (—1/3,1] which means 
that G is a growing time function, G = Go (< + to)" , 
a > 2, e (-1/3,1]. 



2. Perfect fluid solution 

As we already know, the behaviour of the solution is 
given by Eq. ([7]) , so it only remains to know the value of 
the constants (ai)j^2 parameter of the equa- 

tion of state. We have found the following solution for 
the FE (|40ll44p with the conservation equation, div T — 0, 
where the stress-energy tensor is defined by Eq. ([5]) , 

oi = a2 = a3 = 1, q = 0, w = -i, (49) 



C. Scalar models 

1. Scalar model 

The FE for this model are described by Eqs. (1401144^ 
with the stress-energy tensor is defined by Eqs. (I12j|13p 
and the conservation equation given by Eq. (jl4p . Hence 
taking into account Eq. (|15p we find the following result 

ai 02 = a3 = 1, q^O, 
Vme (-1,1)\{0}, a = /? = 2(l-TO2). (51) 
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With these results the metric collapses to Eq. p8|) . so it 
is not inflationary and isotropic. As it has been proven in 
(see for example [l^ and for a extensive review of results 
[7] and all the references therein) we cannot say anything 
about its dynamical behaviour since the obtained value 
of a scape to this study. Nevertheless the study of the 
solution through perturbations shows us that it is stable. 
The potential behaves like the dynamical cosmological 
constant V ^ k ^ t^^. 



2. Scalar model with G-var 

In this case, the model is described by the FE (|40t - 
|44| with the stress-energy tensor defined by Eq. ()12t - 
I13p while the conservation equations is defined by Eq. 
p6)) . The solution behaves as ([17]). We find the following 
solution 



0,2 = as = 1, q = 0, 



(52) 



Vto G (— 1, 1) \ {0} , therefore the metric collapses to Eq. 
(|48| . so the solution is isotropic and not infiationary, 
while 



a 



Go 



2(1- 



> 0. 



(53) 



2a 



With these results we do not know how G ^ (t + <o) 
may behave. Note that if a g (—1, oo)\ {0} , the poten- 
tial, Vit) ^{t + to)"^^"+^^ , behaves as a decreasing time 
function, as a positive dynamical cosmological constant, 
but G is decreasing if a e (—1,0) , and it behaves as a 
growing time function ifi^ a > 0. 



3. Non-interacting scalar and matter model 

The geometric part of the FE for this model are given 
by Eqs. (|40M4|) with the stress-energy tensor defined by 
Eqs. ([6] and [T2lfT3|) while the conservation equations are 
defined by Eqs. (divT™ = = divT"^). The solution 
behaves as (fT8|) . We find the following solution 



(54) 



ai = 02 = as = 1, (7 = 0, 
a = 13, po = 3n-m^-^ 



as it is observed the metric collapses to Eq. (|48| . so 
the solution is not inflationary but it is isotropic and 
only valid for w = — i. We have found the restriction, 
/3 < 2 (l - m2) , Vm e (-1, 1) \ {0} . Therefore a = f3 e 
(0,2). 



4- Non-interacting scalar and matter fields with G-var 

The FE for the model are given by Eqs. (|40ll44p . The 
stress-energy tensor is defined by Eqs. ([5] and [T^]fT5)) 



while the conservation equations are defined by Eqs ([20l 
[^Tjl . The solution behaves as (|22p . therefore we have 
obtained the following solution 

ai=a2 = a3 = l, q = 0, Vm G (-1, 1) \ {0} , 



so the metric collapses to Eq. (|48)) . i.e. the solution is not 
infiationary and isotropic and it is only valid for a; = — . 
The rest of the parameters are 

3 



13 = a' 



Po 



Go 



(1 - TO^) - 3a^ 



(55) 



so, in order to get po > 0, we have the following restric- 
tion, < ■ Unfortunately with these results we 

do not know how G ~ (< -I- ^o)^" may behave. Only by 
"assuming" that the potential, V{t) ~ (i + io)^^^"^^\ 
mimics the dynamic of A we may bound the value of 
a, that is, a € (— l,oo)\{0}, and therefore, as above, 
we may say that G, is decreasing if a £ (—1,0), and it 
behaves as a growing time function iff a > 0. 



5. Interacting scalar and matter model 

a. Approach 1. The FE are given by Eq. (|40ll44p 
and the stress-energy tensor defined by Eqs. ([6l and [121^ 
[T5)) The conservation equation by Eq. (P^ . Since the 
behaviour of the quantities are given by Eq. (jlSp . then 
we have found the following solution 

fli = 02 = as = 1, g = 0, 



a 



2 (1 - w?) - Po (c^ + 1) , 



/3 = 2(l-m2) + ^(c.-l), 



(56) 



Vm e (-1, 1) \ {0} , and Vo; e (-1, 1], with the restric- 

2(1- 

tion Po < ^— 



(cu+l) 



With these results the metric col- 



lapses to Eq. pS)) . therefore the solution does not inflate 
and it is isotropic. For example, if = 0, (which may 
stands for a mixture of dark matter and dark energy) 
then 



/3 



1 



2' Po-2(l 



with a < 2 (l - m^) e (0, 2) . 

b. Approach 2. We have found the following solu- 
tion for the FE (|40j|44p with the stress-energy tensor is 
defined by Eqs. (|6]and fT2llT3l) and the conservation equa- 
tions given by Eqs. (|24ll25p . The behaviour of the quan- 
tities is given by Eq. (IT51) . so 



ai 



Po 



a2 = as 
2 (1 - 
2 (l - rr 



1, q = Q, 

) — a , 

, a<2(l-7 



m ) — a) 



1 



S = ^(Suj + l) 



(57) 
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where to > —1/3, and Vm G (—1, 1) \ {0} , otherwise S > 
(reniember that in this case S must be negative) . With 
these results the metric collapses to Eq. (|48|) so it is 
isotropic and non-inflationary. For example if we set 
w = 0, then the solution reduces to: 



/3 = 1 — rn^ H , po = 2 (l — m^) — a, S = 

2 Vq 



with a < 2 (l - m^) . 

c. Approach 3. The FE (|40ll44p with the stress- 
energy tensor is defined by Eqs. ([5] and [T^lIT^ while 
the conservation equations are given by Eas. (|26ll27|) . As 
in the above studied cases the behaviour of the quantities 
is given by Eq. ([T8)). so we get the following solution 



fll 


= 02 = 




2(1- 






p 


= 2(1- 


5 


1 






^3+' 



0, 



bJ + l 



< 2 (1 - m^) 



(1-^) 



(58) 



and therefore the solution is isotropic and non- 
inflationary since the metric collapses to Eq. (^5)) . As 
it is observed, we have obtained the same solution as in 
the approach 2. In this model, there is a critical value for 
w, w = -1/3, thus uj > -1/3, and Vm e (-1,1)\{0}, 
otherwise 6 < (remember that in this case i5 must be 
positive). As in the above studies cases if we set w = 
then we get 



with Vm G (— 1,1)\{0}, so the solution is non- 
inflationary and isotropic. The rest of parameters are 



1 



GqPq 
1 



[2 (1 - m^) - Go (a2 + po)] 



/3 = ^ [6 (1 - "i') - Go («2 + 2po)] , 
[6 (1 - m^) - Go (3^2 + 2po)] 



3GoPo 



(59) 



therefore the solution is valid for Vui E (—1, 1]. 
If we set w = 0, then we obtain: 



po = 7^ [2 (1 - m^) - Goa^] , 



^ [2 (1 - m^) + Goa^] , 5=' 



3' 



2 ^ 2(l-m^) 



since 



so we have the following restriction, a < — ^ 
po > 0. With these results we do not know liow G ^ 
{t + to)^°' may behave. Note that if a S (— 1, cx))\ {0} , 

the potential, V{t) ~ (i -I- to)~^^"^"^' , behaves as a de- 
creasing time function, as a positive dynamical cosmo- 
logical constant, but G is decreasing if a € (—1, 0) , and 
it behaves as a growing time function iff a > 0. 



D. Scalar tensor model 

In this case the effective stress-energy tensor takes the 
following form 



/3 = 2 (l - m^) + 



a 
2' 



Po 



2 1 



1 

3' 



with a < 2 (l - m^) . 

As it is observed for uj — we have obtained the same 
results in the three approaches. The main difference with 
regard to the non-interacting case is that in this case 
the solution is valid for Vo; e (—1/3,1], except in the 
approach 1 where lo G (—1, 1]. 



6. Interacting scalar and matter fields with G-var 



Tl 



Tl 



-P 



Svr 0' fd' b'\ UJ f(j)' 



d' 



' (j)^ (I) \a b ) 2'^ 



LO (<i>' 

2 



(j) \d a 
Stt 0' fa' b'\ uj f <j)' 



(60) 

^.A(0) 
(61) 

-^ + A(^) 
<P 

(62) 

^.A(0) 
<P 

(63) 



The FE for this model are given by Eqs. (|40j|44p . The 
stress-energy tensor is defined by Eqs. ([S] and [T^lfT^ 
while the conservation equations are given by Eqs. (l29l - 
l30l) . In this case the behaviour of the quantities are given 
by Eq. (|32l) . so we get the following solution with Q = 



and conservation equations 
i3 + 2.mi^-+H^]-2(A- '''^ 



d(t) 



= (P - 3p) , 
(64) 

p' + ip + p)H^O. (65) 



ai = 02 = 03 = 1, q = 0, 



For the FE (|34II35I) and taking into account Eq. 
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we find the next solution. 

ai = a2 = 03 = 1, (/)o = 1, m = m, n — —1 — 87, 
po = -^^(^ ((- + 1) (37 + 1)^ + 67 + 2m^) , 

(((7 - 1) (37 + 1)') ^ + 2 (37 + 1) K - 1)) 



Ao 



2(1 + 7) 



(66) 



therefore, the metric collapses to Eq. (|48|) , so the solution 
is isotropic an non inflationary, po = 0, iff 



7 = 



— m^) + 3 — 2™^ + 2 



3(a; + f) 



7 = , W = 



while, fi = jc ~ If 7 < ^1/3 then n > and 
n<0V7G (-1/3,1]. 

For example, if we set uj = 3300 ([22]) and m = 
±1/2, then po > <^=> 7 G /i, where Ii — 
(-0.34054, -0.32633) , but if m = ±50, then po < 0, 
V7 e (-1,1], and if m = ±0.001, po > 7 G /q, 
with lo = (-0.34164,-0.32523). While Aq > 0, V7 G 
(-l,l]\/2, where h = (-0.333446,-0.333333). As it is 
observed Aq = 0, if 
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1 

'3' ^ 



— f cj ± J2u (2w - 3m2 + 3) 
3a; V 



note that Aq < 0, V7 G /2. 

Therefore this solution has only physical meaning V7 G 
/i, i.e. a small neighborhood of 7c, Si'fc = — f) = 
note that I2 d h- For example for 7c, then po > but 
Aq = and n = 0, therefore (j) = 4'a ^-nd this means 
that G = Gq (i.e. constant), i.e. for this specific value, 
the solution is the same than the one obtained in the 
classical model for a perfect fiuid. In the same way we 
may compare the JBD solution with the one obtained 
in the classical situation when G and A vary. If 7 > 7c, 
V7 G h then h <0 and therefore GeS ~ 4>~^,isa growing 
time function, and A behaves as a positive decreasing 
time function. Note that n > 0, V7 G (-0.34164,7c), 
n = if 7 = 7c, and n < V7 G (7c, -0.32633). 



the following metric 

ds"^ ^ -dt"^ + a^dx^ + (6^ cos^ x + d^ sin^ x) dy'^+ 

2 cos a; sin x{b'^ — d^)dydz + {b^ s\v? x + d^ cos^ x) dz^. 

(68) 

The homothetic vector field is obtained from Eq. ([1]) 

is 

n = {t + io) dt + {l- 02) ydy + (1 - 02) zd,, (69) 

so the scale factors must behave as 

a^ao{t + to), b = bait + taf' , d^do{t + t^T , 

(70) 

with ai = 1,02 = 03 G M^. We emphasize the fact that 
we have been able to obtain non-singular scale factors. 
So the restrictions are: 

ai = 1, 02 = 03, (71) 

hence the metric collapses to the following one 

ds^ = -dt^ + {t + tof dx^ + {t + t^f^ {dy^ + dz^) , 

(72) 

thus it belongs to a localy rotationally symmetric (LRS) 
BVIIo which looks like a LRS BI. In this way we find a 
new KVF, iA=d^. 
The FE read: 

a'b' a'd' b'd' 1 / d^ b^\ 

(73) 

b" d" b'd' I f ^ „ ^^2 

b d b d Aa^ \ b^ d^ ^ 



(74) 

= 0, (75) 



2 + ^-^^- —3 
a b ' a b ' Aa"^ \ d"^ 6^ 



6 dab ad a'^ \ d^ 6^ , 

+ — |2 + — - — ) - -SttGT^, 

(76) 

a" d" a'd' If d^ 3b^\ ^ 4 

a d a d Aa^ \ b^ d^ J 

(77) 

and the conservation equation(s) for the different studied 
cases. 



IV. BIANCHI VIIo 
A. The metric 

We have the following Killing vector fields ([21]): 

= dy, ^2 = 9^, ^3 = 9^- zdy + yd:,, (67) 

with [^1,^2] = 0, [^2,6] = -Ci, = 6, so we may 

consider that Cf^ = 1 , C23 = — 1 , and therefore we have 



B. Classical solutions 

1. Vacuum solution 

There is a solution with 02 = 0. Therefore the metric 
collapses to this one 

ds'^ = -dt'^ + {t + tof dx^ + dy^ + dz^ . (78) 

This metric is the Taub form of a flat space-time (see 
chaper 9 Eq. (9.6) of [H and 
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2. Perfect fluid solution 

The FE are given by Eqs. (|73ll77p with the conser- 
vation equation, divT = 0. The stress-energy tensor is 
defined by Eq. ([6]), and the general behaviour for the 
solution is given by Eq. ([7]). We find only two solutions 

1. This solution is only valid if lu = I, 

ds^ = -dt^ + {t + taf dx^ + dy^ + dz^, (79) 

which is the same metric as the obtained one in the 
vacuum solution. 



2. 



ds^ = -dt^ + {t + tof {dx^ + dy^ + dz^) , (80) 

i.e. ai = 02 = 03 = 1, g 0, a flat FRW and only 
valid for w = — i . A simple calculation shows us the 
obtained solution is non-inflationary, (7 = 0, and 
isotropic since the anisotropic parameters {A, W^) 
vanish. 



3. Time-varying constants scene 

In this case the behaviour of the solution is given by Eq. 
([TT|). The FE (|73ll77p with the conservation equations 
([T0|. where a — {uj + l)h and h = {1 + 2a2). and the 
stress-energy tensor is defined by Eq. Therefore we 
have found the next solution by taking into account the 
obtained SS restrictions for the scale factors given by eq. 
dTlT) . G behaves as follows: 



G — Gq {t + to)"' ^ , Go 



A 



Anpo (w -f 1) a ' 



(81) 



where A = 2a2 + 02 , while the cosmological "constant" 
behaves as: 



A^A I 



^ \t + to) ^ = Ao{t + to) ^ 
a ' 



We find the following solution 

ai = 02 = 03 = 1, g = 0, p = po{t + ta) 
G = Go{t + tof''+'^-\ A = Ao{t + toy 



-3(w+l) 



(82) 



valid Voj £ (—1, 1], and therefore the metric collapses to 
Eq. (|M|) . We find the following behaviour for G and A 

decreasing Vw e [—1,-1/3) 
G ~ ■{ constant if w = — 1 /3 
growing Vw G (—1/3, 1] 

negative Vw e [-1,-1/3) 
An « ^ vanish if = — 1 /3 

positive Vw e (—1/3, 1] 



As it is observed the solution is valid Vw € (—1,1], instead 
of only for w = —1/3. For this value of EoS the solution 
collapses to the PF solution, i.e. G — const, and A = 0. 
The metric collapses to the given one by Eq. ([50)) so it is 
non-inflationary and isotropic. As in the case of the BV 
solution, only under the assumption of A > ([l^), we 
may say that G behaves as a growing time function. 



Scalar models 



Scalar model 



The FE are given by Eqs. (j73ll77|) while the stress- 
energy tensor is defined by Eqs. (|12m3|) . For this model 
the conservation equation is given by Eq. (ITi)) . Hence, 
by taking into account Eq. ((T5|) we find the following 
results 

1. Vacuum solution 

ai = 1, 02 = 03 = 0, a — (3 ^ 0, (83) 

therefore we have obtained the same solution as in 
the vacuum case. 

2. A non-trivial solution 

01=02=03 = 1,17 = 0, a = /3 = 2. (84) 

therefore for this solution the metric collapses to 
Eq. ([M]) i.e. the solution is non-infiationary and 
isotropic. From the DS point of view the solution 
is stable ([7]). 

2. Scalar model with G-var 

The FE of this model are given by Eqs. (|73ll77p with 
the stress-energy tensor is defined by Eq. (|12j|13l) while 
the following conservation equations is defined by Eq. 
(IT51) . The solution behaves as 

= 00 (i + toV" , V{t) ^p{t + to)"'^"+'^ , 
G = Go (t + to)'" ■ 
We find the following solution 



oi = 02 = 03 = 1, g = 0, 
a = l, /3 = 1, Go = 2. 



(85) 



As in the previous studied models, the metric collapses 
to Eq. (IMl) . so it is non-inflationary and isotropic. As 
it is observed, since a = 1, then the potential V is pos- 
itive and behaves a.s V ^ t^^, while G ^ t^, i.e. G is a 
growing time function. Note that GV ^ t~^, while in the 
standard approach A ~ t~^. Note that in the classical ap- 
proach A is not multiplied by G while in the quintessence 
approach the potential V, is multiplied by G. 
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Non-interacting scalar and matter model 



The model is described by the following FE (|73ll77p 
with the stress-energy tensor is defined by Eqs. ([51 and [T^ 
I13p while the following conservation equations is defined 
by Eqs. (divT™ = = divT'^). The solution behaves as 
P^ . We find the following solution 



ai = 02 = as = 1, 9 = 0, 
a = /3 = 2(l-^), poe(0,3), Lo 



(86) 



With these results the metric collapses to Eq. ((5(I| and 
therefore it is non-inflationary and isotropic. The so- 
lution is only valid for u) = — i, with the restriction 
Po G (0, 3) so if po ^ 1 then the solution collapses to 
the one obtained in the scalar model with a = /3 = 2. 



4- Non-interacting scalar and matter fields with G-var 



The FE are given by Eqs. (j73ll77p with the stress- 
energy tensor defined by Eqs. (jH and [T2llT3|) and the 
following conservation equations given by Eqs (|20II21|) . 
The solution behaves as ((^^ , therefore we have obtained 
the following solution 



ai 



0-2 ^ as ^ 
1=13, Po 



1, q 
3 



- 0, 
3 
2' 



(87) 



In order to get, po > 0, we have the following restriction: 
Go < 2. The metric collapses to Eq. ([50)1 and therefore it 
is non-inflationary and isotropic and only valid for uj — 
— i. As in the scalar solution with G— var, since a = 1, 
then the potential V is positive and behaves as ^ ~ t~'^, 
while G ^ , i.e. G is growing time function. 



Interacting scalar and matter model 



b. Approach 2 We obtain the following solution for 
the FE (|73ll77p with stress-energy tensor is defined by 
Eqs. ([6l and [T2lfT3| and the following conservation equa- 
tions given by Eqs. (|24ll25p . The behaviour of the quan- 
tities are given by Eq. ([18]), 



oi = a2 = 03 = 1, g = 0, 



a = 2-po(c^ + l), /? = 2+-po(w-l), 
l + 3a; 



y/2 -po{Lo + iy 



(89) 



valid Vw € (—1/3, 1], (5 < 0) , with the same restriction, 
Po < -(Z^T)- "^^^ metric collapses to Eq. ((5(1)) so it is 
isotropic and non-inflationary. For example if we set 
a; = 0, then we get 

1 a 
<5 = 1=, P^l + TT' Po = 2-q;, 

Va G (0,2). Observe that we get the same result as in 
the previous case. 

c. Approach 3 In this case the FE are given by Eqs. 
(I73II77P while the stress-energy tensor is defined by Eqs. 
(l6land fT2llT3| and the conservation equations are given by 
Eas. (|26ll27)) . As in the above studied cases the behaviour 
of the quantities are given by Eq. (|18p . so we get the 
following solution 

oi = 02 = 03 = 1,9 = 0, 



a = 2 - Po (t^ + 1) , 6 = -+UJ, 
/5 - 2 + ipo (^ - 1) , 



(90) 



valid Vlo £ (—1/3,1], {6 > 0)with the restriction, po < 
j^j^TTj- The metric collapses to Eq. (IM|) and therefore it 
is non-inflationary and isotropic. For example if we set 
id = 0, then we get 



a. Approach 1 The model is described by Eqs. (j73I 
[77)) where the stress-energy tensor is defined by Eqs. ([6] 
and [TWr5)) and the following conservation equation (1^51) . 
Since the behaviour of the quantities are given by Eq. 
(fT5)) . then we have found the following solution 

oi = 02 = 03 = 1, 9 = 0, 

a = 2- poiu + l), (3 = 2+^po{uJ-l), (88) 

this solution is valid Vw G (—1,1], instead of a unique 
value of a; as in the case of non-interacting fluids. We find 
only a restriction, po < (^J^i) ■ The metric collapses to Eq. 
([80]) and therefore it is isotropic and non-inflationary. For 
example if we set ui = 0, then we get 

a = 2(/3-l), po = 2(l-/3), V/3e (1,2). 



1 a 
(5=2, ^ = 1 + -,po = 2 - a, 

Va e (0, 2) . 

As it is observed we have obtained the same solution 
in the three cases. The main difference with the non- 
interacting case is that for these models the solution is 
valid for all values of a;, i.e. Vw G (—1/3, 1], which is more 
realistic. In approach 1 we have obtained lo G (—1, 1]. 



6. Interacting scalar and matter fields with G-var 



The FE (j73ll77p with the stress-energy tensor is defined 
by Eqs. (0 and [T^lfT5|) while the conservation equations 
are given by Eas. (|29ll50)) . In this case the behaviour of the 
quantities are given by Eq. ([5^ . so we get the following 
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solution with Q = SHp„ 



ai = a2 = 03 = 1, 9 = 0, 
1 



GqPq 



[2-Go(l + po)] 



[6-Go(3 + 2po)]. 



3GoPo 



(91) 



The metric collapses to Eq. (|5(71) so it is isotropic and 
non- inflationary. As in the previous cases, since a — 1, 
then the potential V is positive and behaves as V ^ t~'^ , 
while G ^ t^, i.e. G is growing time function. If we set 
w = 0, then we obtain: 

a = l, Po = 7^[2-Go], 

in order to get, po > 0, we have the following restriction 
Go <2. 



D. Scalar tensor model 

The effective stress-energy tensor takes the following 
form 



b' uj ((j)' 



A(^), 



d h 

^1 = — TV 



(92) 
(93) 



^3 =-— P 



(94) 

+ A(0) 
(95) 

A(0) 
(96) 



" 7^ "6/^2 



d a / 2 



ai = 02 = 03 ==1, 9 = 0, 
00 = 1, n = -1 - 37, 



Po 



Ao = 



w (37+ 1)^972 + 127 + 1 



87r(l + 7) 



0. (1-7) (37+1)^ +2(37 + 1) 



2(1 + 7) 



(99) 



With these results, the metric collapses to Eq. (1801) and 
therefore it is isotropic and non-inflationary. Therefore, 
Po = 0, iff 



(w ± V2Z7+3 + 2) , , 1 

7=-^ n — -T, ^^7^-1, ^7^7:> 



3a; + 3 



while, Aq = 0, iff 



7 = -3,7 = 



Lu ± ^2uj (2w + 3) 



'duj 



If we fix w = 3300, ([22]) then we get 



po > 
Ao < 



76/1, /i (-0.3416,-0.3252), 
76/2, /2 = (-0.3334, -0.3333), 



in fact Ao > if 7 ^ /2. Therefore this solution is only 
vahd V7 e /i, a small neighborhood of 7c, f (7c = —\) = 
Ii , note that I2 C ^ijbut if we take into account the cur- 
rent observations (19]) which suggest us that the cos- 
mological constant must be possitive then we arrive to 
the conclusion that the solution is valid V7 G /i\/2- If 
7 > 7c, V7 G /i then n < and therefore Goff « (j)^^ , is 
a growing time function, and A behaves as a positive de- 
creasing time function. For example if 7 = —1/3, which 
means that n = 0, then 



Po 



Ao = 0, Gcff = Go, 



compare with the perfect fluid solution and with classical 
model with TVC. Note that n > 0, V7 G (-0.3416,7c), 
n = if 7 = 7c, and n < 0, V7 G (7c, -0.3252). 



BIANCHI IX 



while the conservation equations are: 
(3 + 2a;) + iJ-^ - 2 A - 0— = — (p - 3p) , 



p' + {p + p)H = 0. 



(97) 
(98) 



For the FE p4ll35l) and taking into account Eq. ((57)) 
we find the next solution 



A. The metric 

The Killings are ((U): 

Ci = dy, 

^2 = cos ydx — cot X sin ydy 
^3 = — sin ydx — cot x cos ydy 



sm y , 



cosy 
sin a; 
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therefore, [^1,6] = 6, [6,6] = -6, [6,6] = 6- Alge- 
bra brings us to get 

ds'^ = -dt'^ + (a^it) sin^ z + b^{t) cos^ z) dx^+ 

2 (6^(t) — a^(i)) sin z sin X cos zdxdy+ 

(a^ {t) sir? X cos^ z + 6^ (i) gin^ a; sin^ z + 

+d^{t) cos^ x) dy^ + 2d'^{t) cosxdydz + d'^{t)d-i^. 

(100) 

The metric admits the fohowing HVF 

n = {t + to) dt + {l- 02) ydy + (1 - as) zd,, (101) 
so the scale factors must behave as 

a^ao{t + to), b^bo{t + to)"' , d = do{t + to)"' , 

and the restrictions on the constants (a,;)^^j^ G M, are 

ai = 1, 02,03 eR. (102) 

For this model the FE read: 

a' h' a! d! V d' 1/1 1 1 
^ 77 i 7 2 W^52 + d2 







62 




^b'^d^ ^ 


02^2 




a" 


6" 




a 


-T^ 






1 






-6^ + 


6" 




6'd' 


T 




^ 77 ' 


U 




d2 




a2d2 + 


62 a2 



(103) 



7 



(104) 



GT4 



(105) 



B. Classical solutions 

1. Vacuum solution 
There is not SS solution for the vacuum model. 



2. Perfect fluid solution 



with the con- 



The FE are given by Eqs. 
servation equation, T/.^- = 0. The stress-energy tensor is 
defined by Eq. ([6|), and the general behaviour for the 
solution is given by Eq. ([7]). We find the following solu- 
tion 

ai = 02 = flg = 1, q = 0, w -1/3. (109) 
In this way the metric collapses to this one: 
ds'^ = -dt^ + 

{t + tof (dx^ + dy^ + dz^) + 2 (t + t^f cosxdydz 

(110) 

and therefore there is another KVF, = dz- As in the 
previous models this solution is isotropic, {A, W^) vanish 
and it is only valid for the EoS uj — —1/3. The obtained 
solution is stable from the dynamical point of view (0). 

3. Time varying constants framework 

In this framework the FE are given by Eqs. (|103I 
llOSp with the conservation equations ([rU]). where a = 
{uj + 1) h and h = (1 -|- 02 -l- 03), and the stress-energy 
tensor is defined by Eq. In this case the behaviour 
of the solution is given by Eq. ([TTt . Therefore we have 
found the next solution by taking into account the ob- 
tained SS restrictions for the scale factors given by Eq. 



0" 


d" 


a' d' 




1 




Oi = O2 








^2(02- 




52) 





' d 


a d 




^2- 







U- 



15 



+ (100) 



487rpo (w + 1) ' 



A = il(l 



4c2 V 3(w + 1) 



e (-1,1] 



(t + to) \ (111) 



1 



1 



06 ad 

62 d2 



as above the metric collapses to Eq. (|110p . so it is 
isotropic and non-inflationary. The behaviour of the 
"constants" is the following one: 



0, 



62 rf2 ^2^2 52^2 

" 6" o' 6' 1 / 1 1 1 



(107) 



06 2 V o 



62 rf2 



G 



An 



4 V^^d^ o2d2 62a2 
and the conservation equation(s) 



4 , 



(108) 



decreasing Vw e [—1, —1/3) 
constant if w = — 1 /3 
growing Voj e (—1/3, 1] 

negative e [-1,-1/3) 
vanish if w = —1/3 
positive Voj e (—1/3, 1] 



i.e. we have found the same behaviour as in the previous 
studied models. 
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C. Scalar models 

1. Scalar model 

The FE are given by Eqs. (|103m08|) while the stress- 
energy tensor is defined by Eqs. (I12j|13p . For this model 
the conservation equation given by Eq. (jl4[) . Hence, by 
taking into account Eq. (jl5|) we find the following non- 
trivial solution 

5 

ai = a2 = as = 1, « = /?=-, (112) 

therefore the metric collapses to Eq. (jllOp and the solu- 
tion is isotropic and non-inflationary. As it is showed in 
(0) tfie solution is stable from the DS point of view. 



the corresponding conservation equations given by Eqs 
(|20j|2ip . The solution behaves as ([2^ . therefore we have 
obtained the following solution 

oi = 02 = as = 1, q = 0, 

Po = 4^ (5 - 2Go) , ^ = (115) 

so we have the following restriction, Gq < |, in order 
to get pq > 0. The metric collapses to Eq. (|110p so 
the solution is isotropic and non-inflationary and only 
valid for w = — i. Note that a = 1, so G is a growing 
time function, G ^ , while the potential is positive and 
behaves as a decreasing time function. 



2. Scalar model with G-var 



5. Interacting scalar and matter model 



For this model the FE are given by Eqs. (jlOSHlOSI) . 
The stress-energy tensor is defined by Eq. ^21113^ while 
the conservation equation is defined by Eq. ()16p . The 
solution behaves as ([T7)) . We find the following solution 



ai 



^2 



as = 1, 



q 



a=l, (3 = 1, Go = 



(113) 



where, as it is observed, the metric collapses to Eq. 
(jllOp and therefore the solution is isotropic and non- 
inflationary. Since, a = 1, then we may say that G ~ t^, 
i.e. it is a growing time function while the potential is 
positive and decreasing, it behaves as, V ~ . 



3. Non-interacting scalar and matter model 

In this case, the FE for the model are described by 
Eqs. (|103m08p . and the stress-energy tensor is defined by 
Eqs. (|6l and [T2]fT3l) while the conservation equations are 

= (T"^)^.,). The solution 



defined by Eqs. ((T" 



behaves as (llSp . We find the following solution 



ai = a2 



as 



Po 



1, q 
3/5 
~ 2 V2 ~ 



0, 



(114) 



therefore the metric collapses to Eq. (jllOp and the solu- 
tion is non-inflationary and isotropic. As in the previous 
models, the solution is only valid for cj = 
parameters (a, /3) are only bounded. 



i, while the 



4- Non-tnteractmg scalar and matter fields with G-var 



The model is described by the FE 



with the 



stress-energy tensor deflned by Eqs. (0 and [l^llT^l) and 



Approach 1 The model is described by the FE 
l llOSp with the stress-energy tensor is defined by Eqs. 
(0 and [T^lfT^ and conservation equation (^5]) . Since the 
behaviour of the quantities are given by Eq. (jlSp . then 
we have found the following solution 



ai = a2 = as = 1, 
a = o - Po (w + 1) , 



9 = 0, 

P=\ + \p^{^-l), (116) 



finding that the solution is valid Vcj g (—1,1], with the re- 



striction, Pq < 



.The metric collapses to Eq. (|110p 



2(w+l)- 

which means that the solution is non-infiationary and 
isotropic. For example if we set w = 0, then we get 



a = 2/3- 



Po 



2(3, y(3e 



b. Approach 2 The FE for this model are described 
by Eqs. (|103m08|) with the stress-energy tensor is defined 
by Eqs. (|6land [T2]fT3)) and the conservation equations are 
given by Eqs. (j24j|25p . The behaviour of the quantities 
is given by Eq. (IT51) . so we get the following solution 

ai = 02 = as = 1, <7 = 0, 
5 , „ 5 



a 



2 - po (w + 1) , /3 = 2 + 2/^0 (w - 1) 



6 = -- 



2(1 + 3w) 



(117) 

2{:j+T)' € (-1/3,1], 

((5 < 0) . With these results, the metric collapses to Eq. 
()110p and therefore the solution is isotropic and non- 
inflationary. For example if we set uj = 0, then we get 



Vl0-4po {oj + iy 
finding the restriction, po < ^ 



1 



^ 4 2' 



»4 



Va e 



and therefore we have obtained the same result as in the 
approach 1, as it is expected. 
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c. Approach 3 The FE (|103m08p with the stress- 
energy tensor is defined by Eqs. dH and [T^lfT^ and con- 
servation equations given by Eqs. ()26ll27l) admit the fol- 
lowing solution, note that the behaviour of the quantities 
is given by Eq. ([T8| , 



ai 
a 



as = 1, 



0, 



a2 



5 , s c 1 

2-2(^ + 1), ^=3 



(118) 



with po = Pq, finding the restriction, po < 2{J'+i} ' ^ 
(— 1/3,1],((5>0). Therefore the metric collapses to Eq. 
(jllOp which means that the solution is non-inflationary 
and isotropic. If we set lo — 0, then it is obtained the 
following result: 



Po 



Va e 0, 



and therefore, once again we have shown that the three 
approaches are identical. 



6. Interacting scalar and matter fields with G-var 



D. Scalar tensor model 

The effective stress-energy tensor takes the following 
form 



1 Stt 6' u ( 6'\ , , , 



^2 = — T'P 



(119) 

A(0), 
(120) 



d h 2 



^3 37r 0' (d! , a'\ uj f<j)'Y 4>" 



(121) 



= — rP-^ - + T -ohr -^ + A(<^), 



a h 2 



with the conservation equations 



(122) 



(i + 2u^m[^+H^]-2{k-<^'^ = ^-^{p-Zp), 

(123) 

p' + {p + p)H = Q. (124) 



For the FE (|34ll35p and taking into account Eq. ([32] 
we find the next solution 



The model is described by the FE (|103m08l) with the 
stress-energy tensor defined by Eqs. (|6] and [T2]fT3l) while 

In 



the conservation equations are given by Eqs. (|29ll30l) 
this case the behaviour of the quantities are given by Eq. 
(|32|l . so we get the following solution with, Q = SHpm, 



fli = 02 = 03 = 1 
1 



q — 0, a = 1, 
[5-2Go(l + po)], 



2G0P0 

f3=-^[l5-2Go{l + 2po)], 

S = [15 - 2Go (3 + 2po)] 

6G0P0 



The metric collapses to Eq. (|110p . so the solution is 
isotropic and non-infiationary. Since a — 1, then G is 
growing, it behaves as G ~ i^, while the potential be- 
haves as a positive decreasing time function. Setting 
w = 0, which may be interpreted as a model with in- 
teracting DM and DE, then we obtain: 



a = 1, Po = [5 - 2Go] , 
/3 = ^[5 + 2Go], 6=1, 



so we have the following restriction, Go < |, in order to 
get Po > 0. 



ai = a2 = as = 1, 9 = 

^2 , -,0 9 



Po 



Ao 



2w(37-f 1) 



<f>o = 1 
187^ -h 247 + 1 



167r(l-|-7) 
2u ((1 - 7) (37 + 1)^) + 5 (37 1) 



4(1+7) 



(125) 



therefore, the metric collapses to Eq. (IllOp . so the solu- 
tion is isotropic and non-infiationary. po = 0, iS 



7 



2w ± ^2 (5w + 7) + 4 



6(0; -Hi) 

and the cosmological constant vanish when, Ao = 0, iff 



7 



1 

'3' ^ 



2uj ± y/2uj (8a; 4- 15) 



6u; 



If we fix a; = 3300, {^) then we get 



Po>0 
Ao <0 



76/1 Ii = (-0.3426, -0.3242) , 
7^/2 h = (-0.3335, -0.3333), 



in fact Ao > if 7 ^ /2 . Note that h <Z h - Therefore this 
solution is only valid if 7 G /i , i.e. for a small neighbor- 
hood of 7c, £{^c — —\) — Ii, note that I2 C Ii- As in 
the above cases, if we consider that Aq must be possitive 
then we arrive to the conclusion that the solution is valid 
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V7 € /i\/2- If 7 > 7c, V7 G /i then n < and therefore 
Geff ~ (/)~^, is a growing time function, and A behaves as 
a positive decreasing time function. Setting 7 — —1/3, 
then n = 0, and therefore 



i.e. we obtain the classical perfect fluid solution. 

VI. CONCLUSIONS 

In this paper we have studied how the constants G 
and A may vary in different theoretical models as general 
relativity with a perfect fluid, scalar cosmological mod- 
els ( "quintessence" models) with and without interacting 
scalar and matter fields and a scalar-tensor model with 
a dynamical A. We have applied the outlined program 
to study three different geometries which generalize the 
FRW ones, which are Bianchi V, VIIq and IX, under the 
self-similarity hypothesis. These geometries are in prin- 
ciple homogeneous but anisotropic. We have put special 
emphasis in calculating exact power-law solutions which 
allow us to compare the different models. 

As we have shown, in all the studied cases, we arrive 
to very similar conclusions. The first of them is that the 
solutions are isotropic (^,>V^ vanish) in spite of consider- 
ing anisotropic geometries, and noninflationary ((7 = 0). 

With regard to the classical solutions i.e. solutions ob- 
tained within the general relativity (with a perfect fluid 
as matter model) and its TVC model, we arrive to the 
conclusion that the obtained solutions are only valid for 
a unique value of the EoS, w = —1/3. In the Bianchi 
V model the metric collapses to a metric with 6 KVF, 
therefore, as we already know, the FE generalize the case 
FRW with negative curvature. In the Bianchi VIIq model 
the metric collapses to the flat FRW and for the Bianchi 
IX the obtained metric has only 4 KVF, nevertheless the 
FE generalize the FRW with positive curvature. In the 
case of TVC models, we have arrived to the conclusion 
that in all the cases the solutions are valid Vw € (—1,1]. 
For uj = —1/3 the solutions collapse to the standard so- 
lution with G = const, and A = 0. If we consider the 
current observations which suggest Aq > 0, then we ar- 
rive to the conclusion that the solution is valid only for 
Vo; G (—1/3,1] and that the gravitational "constant", 
G{t) ^ behaves as a growing time function while 

With regard to the scalar cosmological models, we 
conclude that in all the cases, the metrics collapse to 
the ones obtained for the case of a perfect fluid model. 
All the obtained solutions are stable from the dynami- 
cal system point of view and therefore relevant from the 
physical point of view. In the case of the scalar models 
with G— var, we arrive to the conclusion that in the cases 
BVIIq and BIX, G behaves as a growing time function, 
G{t) ~ i^, and the potential, which mimics a dynamical 
cosmological constant behaves as a positive decreasing 



function, V ^ t~'^. In the BV model, since it depends 
on more parameters like m, then the obtained solution is 
not so precise in such a way that we only may say that if 
a e (-1, oo)\ {0} , the potential, V{t) ^ {t + ^o)"^^"'^^^ , 
behaves as a decreasing time function, as a positive dy- 
namical cosmological constant, but G is decreasing if 
a G (—1,0), and it behaves as a growing time func- 
tion iff a > 0. In the non-interacting case, we arrive 
to the conclusion that the obtained solutions are only 
valid for a unique value of the EoS, ui = —1/3. In the 
non-interacting case with G— var, the solutions are only 
valid for a; = —1/3, while G and V behave like in the case 
of the scalar model with G— var, so we have no new in- 
formation about their behavior. For the interesting cases 
of interacting scalar and matter fields we arrive to the 
conclusion that for all the cases the solutions are very 
similar. For the approaches 2 and 3 the solutions are 
the same and only difference with respect to approach 1 
is the range of validity for cj, the parameter of the EoS. 
The main difference with respect to the non-interacting 
case is that the solutions are valid Vw € (—1/3, 1] in ap- 
proaches 2 and 3 while in approach 1 Voj e (—1, 1]. This 
situation is more realistic from the physical point of view. 
For the critical value uj = —1/3 the solutions collapse to 
the non-interacting case in approaches 2 and 3. For the 
case of interacting scalar and matter models with G— var 
we conclude that the solutions are valid Vw € (—1/3, 1], 
while G behaves as a growing time function G{t) ^ i^, 
and the potential behaves as a positive decreasing func- 
tion, V ^ i""*, as in the above studied models. 

As we have commented along the paper, the SSS may 
be considered as asymptotic solutions for general so- 
lutions. In the same way we have pointed out that 
the PF and the scalar solutions are stable from the 
DS point of view ([2l|-0). Since in the TVC scheme 
we have obtained the same behaviour for the scale fac- 
tors than in these solutions, we may ^'conjecture^^ that 
the solutions obtained in this framework are also sta- 
ble, note that always we obtain relationships like that, 
Gpm ^ t^'^ ^ Gpcj,. A simple perturbation analysis may 
be carried out in order to prove this conjecture. 

For the scalar-tensor theories we arrive in all the cases 
that the solution is only valid for a small a small neigh- 
borhood of 7c, £["ic = — -j) = h- For 7c, then po > 
but Aq = and n = 0, therefore (f) ~ (pQ and this means 
that G — Go (i.e. constant), i.e. for this specific value, 
the solution is the same than the one obtained in the 
classical model for a perfect fluid. In the same way we 
may compare the JBD solution with the one obtained in 
the classical situation when G and A vary. If 7 > 7c, 
V7 G Ii then n < and therefore Goff » is a grow- 
ing time function, and A behaves as a positive decreasing 
time function. Note that n > 0, V7 G f , 7 < 7c, n = if 
7 = 7c, and n < V7 G £, 7 > 7c. Once again, if we take 
into account the current observations, then we arrive to 
the conclusion that G behaves as a growing time function 
G{t) ~ t^'^~^, and the cosmological constant behaves as 
a positive decreasing function, A ~ as in the above 
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studied models. We have performed the same calculation 
considering the scalar-tensor theory defined by 



1 

8^ 



C 



M 



arriving to the same solutions that the obtained ones 
where, which is not obvious as we have pointed out in 
(i). 
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Appendix A: Interacting scalar and matter models 

In this appendix we will prove the results stated in 
section II. We shall study the admitted form for the po- 
tential, as well as, for the energy density, in the SS frame- 
work. For this purpose we study the different DE through 
the Lie group method ([Til). 



1. Approach 1 



(Al) 



where r = {cu + 1) h, H = ht ^. Lie's method applied to 
Eq. (|A1[) brings us to get 



We may rewrite Eq. (1231) in an alternative way 



t £,<t><t> = 0, 

2t'^r]t4, - t^^tt + 3i2^0V0 + th^t - ^ 0, 
t^TlV^xt, + t^Vtt + 4i^0 [tp' + rp) + 
t^V^ (2^4 - ri^) + htT]t = 0, 
tp' {3t^t - 2i?70 + rO + 
rp(3t6-2i??0-O+^V = O, 
T]t {rp + tp') = 0. 



(A2) 
(A3) 
(A4) 



(A5) 



(A6) 



(A7) 



Then, we may impose some symmetries in such a 
way that condition on the integrability of the functions 
V and p are obtained. For example, the symmetry 



= t,?7 = 1) , which brings us to get the invariant so- 
lution (j> — Int, we get, from Eq. (|A5p 



V^^ + 2V^ = 0, 
and from Eq. (|A6[) 



p" = -(3 + r)^-2r^ 



V = exp (-20) , 



p = Cit-^ + C2t-'\ (A8) 



setting 6*2 = 0, we get 

= lni, 1/ = exp (-2(/)) , p^C\r^. 



(A9) 



2. Approach 2 

In this model the conservation equations read: 

p'^ + {lj + 1)p„,H = -cf>'qt (AlO) 

(All) 



If we assume such functional relationship, q'^ = f3pm, 
then we obtain for Eq. (|A11[) the following system of 
PDE 

<2^^^ = 0, (A12) 
-2t2^t0 + 2ht^^ + t^T]^^ = 0, (A13) 
3i'C0 {V^ - I3p) + thit -K + 2t\^ - t^itt = 0, (A14) 



^^77^00 + t^mt + 21^ it (V^0 - /3p) - 



iy^ - Pp) + htr]t - pt'ip' = 



(A15) 



The symmetry (^ = t, 77 = 1) =^ = ln<, brings us to 
obtain the following restriction, from Eq. (jA15[) we get 

1/00 + 21/0-/3 {tp' + 2p) , 

V = exp(-2(/)) = /3p»t-^, 

While Eq. (jAlOp admits the following solutions 



= --(rlni-lnp) + Ci, 



in an alternative way 



p = exp [— — r In t + /3Ci] ~ po 



1 



r + /3 = 2. Therefore we have obtained: 

= lnt, y = exp(-20), p = pot-^. (AI6) 
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3. Approach 3 

The conservation equations read 

p'„ + {UJ +1) PmH = 5Hpm, 

0' (u^+'^\ ^-6Hpm 



(A17) 
(A18) 



From Eq. (|A17I) we get 



and therefore Eq. (|A18[) yields 

dV 



0' ( 0" + h<i)'t-^ + 



-At", 



(A19) 



with a = h{5 — {uj + — A = 6hpo. We have the next 
system of PDE for Eq. (A19\ 



(,4,4,t'^ = 0, 

-2fit4, + 2/lt^0 + t^jy^^ = 0, 

3t^C</."^0 + -K + 2t^Tjt4. - t^^tt = 0, 
3Af+''£,t - 2Af+''r]^ + aAt^+''^ = 0, 
At^+^rjt = 0. 



t^vVc^c^+t\t+^At^+''UV^- 



t^V4.V4, + htrit + 2t^itV4, = 0, 



(A20) 
(A21) 
(A22) 
(A23) 
(A24) 



(A25) 



The symmetry — 1,11 — 1) , bring us to get from Eq. 
(IA25I1 

e'nV^^ + 2e^tVc^ = Q, 

3At2+"^t + aAt^+''^ = 0, 

simphfying 

¥4,4, + 2aV^ = = exp(-2a0), 

where a = — 3. Therefore we get: 

Pm-po^"^ (j^^lnt, F = exp (-2</)) . (A26) 

4. A generalization 

In this model we try to generalize the previous ap- 
proaches. We consider a function Q = 6Hf{t), f{t) is an 
unknown function. Therefore applying the same proce- 
dure we would like to determine the admitted form for 
/ in order to obtain a SSS. For this purpose we consider 
the following system of DE 



p'„, + {Lo + l)p„,H = 5Hf{t), 
0' (u4>+'^\^-6Hf{t), 



From Eq. (|A27|) we get 



p=\5h f{ty-'dt + Ci t 



(A27) 
(A28) 



(A29) 



where r ~ h(uj + \) , and H = ht ^ . In the same way we 
may analyze Eq. (|A28p through the Lie group method 

0' (^0" H'r^ + ^) " -Sht-^f{t), (A30) 

getting the following system of PDE: 

^00^2=0, (A31) 

-2t^^t4, + 2hti^ + t^r,^^ = 0, (A32) 

M^^V^ + th^t -h^ + 2t^r]t4, - t^Ctt = 0, (A33) 

3Shtf{t)^t - 2Shtf{t)f]4, + 6h (tf'{t) - /) C = 0, (A34) 

Shtf{t)rjt ^ 0, (A35) 



fvV4,4, + t^Vtt + 4Shtf{t)C4,V4,- 



(A36) 



The symmetry = t.rj = 1) , bring us to get from Eq. 
(|A36l) 

T^00 + 2F^ = V = Vo exp(-20), 

tf\t) + 2/ = / - /o^-^ 



then Eq. (|A29|) gives 



(A37) 



and therefore 



p = + C2r^ 

Therefore, setting C2 = 0, we have obtained the fol- 
lowing results: 

Pm = Pot~'^, (t) = lnt, V ^exp{~2(j)), / = /ot"^ 

(A38) 

As we can see we can consider any function / that 
behaves as / = fot~^. This fact justifies the employ of 
/ = Pm, as in the previous approaches, but we can also 
consider other forms for / as for example, f — p^ k. t~^, 
or f = Pm + Ptj, ~ t~'^, etc. 
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Appendix B: Matter and scalar fields witli 
G— varying 

In this appendix we will study the equation 



2t^G^ {p + G~-^G'V + t-^rp + G' p) r]t = 0. (B7) 

The symmetry — t, 77 = a0) (j) = t"', bring us 
to obtain the following restriction on the functions G, V 
and p. From (jB4l) we obtain 



G" 



G'2 G' 



G 



Pm + +V , 



(Bl) 



G t ' 
Now from Eq. (IBSp we get 



that we may study as follows. In the first case and in 
analogy with the previous appendix we will study the 
whole equation i.e. Eq. (jBl[) without spliting it while 
in the second case we will split it for a suitable coupling 
function. 



1. Case 1 



Eq. (jBl[) may be written in the following form 



G = Gnta. 



V ((^) = Vof' 



(j}" + ht-' + TTT^ U' + ^ 



2G 



I J 



Qf 



G" 

rt-^ + ) = 0, 



where r = (1 + 7) /i and H = ht ^. The standard proce- 
dure brings us to get the following system of PDE 

t^G'^i^^ = 0, (B2) 

2fG^7j^^ - 4t^G^^t4, + 2iGC0 {2hG + tG') = 0, (B3) In this case we may spht Eq. dHU) as follows 



From Eq. (jB6l) we obtain since cj) — and V (0) = 

Vo(/."'°"" = Vot^^'^-'^\ and G = Ggta then we get 

t 

finding that a solution is: 

p„ = po(t + to)2("-i). (B8) 

Therefore we get (after redefining the constants) 

= 00 (< + to)-'' , Vit) = /3{t + to)-'^"+''> , 
G = Go (< + to)' , - po (t + to)"'^"+'^ , (B9) 

with g = 2a, since Gp w <^^, and p„i = cupm- 
2. Case 2 



etG^e^ (p + tV^,) + 4fG\^ - 2t^G^^t 



-2hG^ - t^G'^ + t^GG") ^ + tG {2hG + tG') = 0, 

(B4) 

ifG^V4,^t - 2t''G^V^r]^+ 

8tG {tG' {V + p)+ rpG + tGp') ^0+ 



2t^G^r]tt + 2t^G^V^^'n + tG {2hG + tG') f]t = 0, (B5) 



G' 

P'm + {P-ni + Pm) H + -q Prn = 6Hpm, (BIO) 

" ' ^'H + + ^p^ = -5Hpm, (Bll) 

in such a way that if we study both equations through 
the LG method then we get. 
For Eq. (IB11|) we get: 

t'G^S,^^ = Q, (B12) 
2t^Gh^^ - 4t^G^£,t^ + 2tG£,^ {2hG + tG') = 0, (B13) 



2t^GG'V4.r] + QtG [tG'V + rpG + tG' p + tGp') 
UG {tG'V + rpG + tG'p + tGp') 1^+ 



2t^G {G"V + G'p' + t-^Grp' + G" p) ^ 



At^G\^ - 2t^G^^tt + 6t^G^V^^^+ 



{-2hG^ - t^G'^ + t^GG") C + tG (2/iG + tG') = 0, 

(B14) 

8iG (tG'V^ + Shp) + 2t^G^r]u + tG (2/iG + tG') rjt 



+ 2t^G{Gp" -rpGt-^ ~G-^G'^iV + p))^ = 0, (B6) + 2t^G^V^^7] - 2t^G^V^r]^ + U^G^V^^t = 0, (B15) 
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2t^GG'V^^r] + 6tG {tG'V^ + G5hp) ft density 
- AtG {tG'V^ + GShp) 770+ 

2{l^a)p + pt = 0, =^ p = pot^'^°''^\ (B18) 

2 (ShG'^ {tp' - p) + 2fV^ {GG" - G'^)) f = 0, (B16) ^^^^ that for Eq. (|BTo1) we get a direct integration since 
2tG {tG'V^ + GShp) rjt = 0. (B17) 

G' 

The symmetry — t, -q = acf)) =^ if) = t"', bring us Pm + (1 + 7 ~ ^) Hpm + 'qP"^ ~ 

to get from Eq. (|5T4)) 

^,2 ^, then pniG = Kt^'', where r = (1 + 7 - 5) . 

G" = -— ==^ G = Got^ . Therefore we have arrived to the same conclusion as in 

G t 

From Eq. (|BT51) 



the above case i.e.: 



</> = <^o + to)"" , V{t) = p{t + to)"'^"+'^ , 
y00 = fc^^ =^ F = G = Go(t + tor, P™ = Po(t + M"''"+'\ (B19) 



Now, from Eq. (jB16p wc obtain an ODE for the energy with g — 2a, since Gp « <^^. 
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